Classical probability waves 

Marius Grigorescu 



Probability waves in tlie configuration space are associated witli coherent so- 
lutions of the classical Liouville or Fokker-Planck equations. Distributions 
localized in the momentum space provide action waves, specified by the 
probability density and the generating function of the Hamilton-Jacobi the- 
ory. It is shown that by introducing a minimum distance in the coordinate 
space, the action distributions aquire the energy dispersion specific to the 
quantum objects. At finite temperature, probability density waves propa- 
gating with the sound velocity are obtained as nonstationary solutions of the 
classical Fokker-Planck equation. The results suggest that in a system of 
quantum Brownian particles, a transition from complex to real probability 
waves could be observed. 
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1 Introduction 



Probability waves, beside classical degrees of freedom, describe the out- 
come of the measurement process in quantum mechanics, or mixed classical- 
quantum systems. However, well beyond atomic scale, quantum coherence 
phenomena appear in macroscopic superconducting devices, while many nu- 
clear collective states are related to classical rotational or vibrational degrees 
of freedom. This wide domain of interplay between classical and quantal in- 
dicates that probability waves could be relevant not only for the atomic 
systems, but also at classical level. 

The early attempts to understand atomic phenomena relied on classical 
analogies, inferring for instance that spectra of hydrogen resemble acoustic 
phenomena in organ pipes [1]. Though, it was the emission, absorbtion and 
thermalization of electromagnetic radiation that provided the interface used 
to shape the general concept of wave-particle duality. 

In the quantum theory, unlike acoustics, the linearity expressed by the 
superposition principle holds only for complex wave functions having a non- 
linear dependence on the observables, such as probability density. Still, the 
underlying classical Hamilton equations and the time-dependent Schrodinger 
equation (TDSE) share a common mathematical framework, as both can be 
obtained from variational principles with suitable action integrals [2j. This 
formal aspect allows to derive the mixed classical-quantum dynamics from 
extended variational principles (HJHllS]. In particular, a quantum particle 
coupled to a classical environment at finite temperature can be described 
using Schrodinger-Langevin [U [6] or quantum Fokker-Planck [5] equations. 
Within this approach, the probability distribution for the quantum compo- 
nent of a mixed system is thermalized by the classical noise, but the quantum 
backreaction has no such effect on the classical component, in the sense that 
no random terms appear in the classical equations of motion as a result of 
purely quantum fluctuations. 

In semiclassical gravity theory, the metric fluctuations appear as ran- 
dom forces with finite correlation time Quantum fluctuations can also 
be included in time-dependent mean-field dynamics, by assigning to the trial 
state an effective temperature related to its energy dispersion, and a "quan- 
tum Langevin force" [8]. Though, such procedures remain ambiguous, as 
the various attempts [9] to understand the transition between classical and 
quantum fluctuations are incomplete and faced with difficulties. 

In this work the probability waves are associated with coherent solutions 
of the classical transport equations, rather than with the continuous action 
of random forces. The property of coherence is attributed to solutions which 
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evolve without changing their functional dependence on coordinates and mo- 
menta. The kinematics of pure distributions, appearing in classical Hamil- 
tonian systems, is discussed in Sect. 2. An important example is provided 
by the action distributions, localized in the momentum space. It is shown 
that by introducing a minimum distance in the coordinate space, the action 
distributions take the form of the Wigner functions. At finite temperature, 
the thermal noise can change the statistical ensemble of a pure distribution. 
This aspect is presented in Sect. 3, considering both stationary and nonsta- 
tionary solutions of the classical Fokker-Planck equation. The results show 
that beside the action and quantum waves, relevant at zero temperature, at 
thermal equilibrium probability waves propagating with the sound velocity 
can appear. Conclusions are summarized in Sect. 4. 

2 Classical coherent states 

The statistical properties of classical systems composed of particles can 
be described by a time-dependent distribution function / > defined on the 
one-particle momentum phase-space M. If 5Qrn is a volume element around 
the point m £ M, then f{m, t)6Qm is proportional to the probability to find 
a particle localized in 5Qrn- Therefore, / is normalized by the integrality 
condition 

J dn^f{m,t) = N , N>1 . (1) 

Let us consider a one-dimensional system consisting of a single particle 
(A^ = 1) with the Hamiltonian H{x,p,t) depending on the coordinate x, 
the canonical momentum p and time. If there are no external forces, the 
equations of motion are 

dtx = dpH , dtp = -d^H , (2) 

where dt = d/dt denotes the total derivative with respect to the time t, and 
dx = d/dX the partial derivative with respect to the variable X. 

The distribution function f{x,p,t) evolves according to the Liouville 
equation 

dtf + CHf = (3) 
where Cnf = ~{H,f} is the Lie derivative defined by the Poisson bracket 

m 

Ch = {dpH)d^ - {d^H)dp . (4) 
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For a Hamiltonian of the form H{x,p) = /2m + V{x), ([3]) becomes 

dtf + -dJ-V'd,f = G , (5) 
m 

where V' = dxV . To find coherent solutions of this equation it is convenient 
to use the Fourier transform f{x,k,t) in the momentum coordinate, 

f{x,k,t) ^ [ dpe'^Pfix,p,t) . (6) 



Thus, if f{x,p,t) is a solution of ([5]) then its Fourier transform f{x,k,t) 
should satisfy 

dtf--dlJ + ikV'f = . (7) 
m 

Various local quantities of interest, such as densities of localization prob- 
ability n{x,t), current j{x,t) and kinetic energy e(x,t), can be expressed 
directly in terms of / and its derivatives at A; = by 

n{x,t)= J dp f{x,p,t) = f{x,0,t) , (8) 

jix,t) = [dp -f{x,p,t) = --dkf{x,0,t) , (9) 
J m m 

e{x,t) = j dp ^f{x,P,t) = -^dlf{x,0,t) . (10) 

In general, f{x,p,t) is specified by the infinite series of partial derivatives 
d^f\k=o, fJ- = 0,1,2,.... Though, certain coherent solution can be defined 
only in terms of n{x,t), or a simple functional of n{x,t). Such functionals 
which satisfy the superposition principle and generate a linear space will be 
refered as " probability waves" . 

An important class of coherent states for the Liouville equation ([5]) are 
the "action distributions" 

fo{x,p,t) =n{x,t)6{p-dxS{x,t)) , (11) 

which remain all the time a product between n{x,t) and 6{p — dxS{x,t)). 
The two real functions of coordinate and time, n{x, t) and S{x, t), are related 
to the Hamiltonian fiow. By Fourier transform (jlip becomes 

fo{x, k, t) = n(x, t)e^^^-^(^'*) (12) 

while ([7]) reduces to the system of equations 

dtn = -dxj (13) 
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na,.[9t5 + %^ + y]=0 (14) 

where j = ndxS/m denotes the current density Q. Thus, by assuming the 
existence of a "momentum potential" S{x,t), we obtain both the continuity 
equation (jl3p and the Hamilton- Jacobi equation [11] in the form (|14p . These 
coupled equations describe "action waves" n{x,t) and can also be derived 
by using a variational approach (Appendix 1). 

The partial derivative dxS{x,t) in (fT2l) is the limit of [S{x + £/2,t) — 
S{x — i/2, t)]/i when f — > 0. If a new parameter a = £/k is introduced, then 

0" J\i 0" k 

Ux,k,t) = Y,raJ*{x-—,t)^{x + —,t) (15) 

where ip{x,t) denotes the complex functional ^ = ■y/nexp(i5/cr)0. Presum- 
ing that, as might be the case in numerical calculations, when both i and k 
decrease to zero the ratio a = i/k remains finite, 

lim 0- > , (16) 

then beside of the limit (|15p we may consider also the "quantum distribu- 
tion" 

/v,(x,A;,t) = r(^- + = (C/-fcV*)(^7fcV') (17) 

as a possible coherent solution of ([7]), presenting its own interest. Here 
i7fc = exp(akdx/2), while o" is a constant. Worth noting, if a = h, then 
obtained by inverting ([6]), 

U{x,p,t) = ^jdk e-'^^~U{x,k,t) (18) 

is the Wigner transform pS] of 'ip{x,t), while the general expressions ([8]), 
dn]) and (jlOp correspond exactly to the canonical quantization of the mo- 
mentum and kinetic energy. Moreover, the probability density v{p, t) of the 
momentum, 

i^{p,t) = J dx f^{x,p,t) = ^ j dx j dk e~''^^f^{x,k,t) , 
can be written in the form i'{p,t) = \ {tpp\'>p)\'^ with 

V'p(x) = ^e^^P/^ , iijpli;) ^ J dx rp{x)i'{x,t) , 



^If 5* is not single- valued, ip is & superposition of terms from different branches [12] , 
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showing that (|17p also contains the interpretation of the scalar product 
(V'pIV') in the Hilbert space TL generated by V' as a probability amplitude. 

The transition from the action distribution /o to the quantum distribu- 
tion could be related to the existence of a limit speed c and a minimum 
distance (. ~ m~^. Thus, if the limit speed is taken into account by restrict- 
ing the integral ([6]) over p to the finite interval [— mc, mc], then (jlSp becomes 
a Fourier series. This series is a sum over an infinite set of discrete values 
of k separated by 1/mc, and therefore a ~ mcH. is finite. 

In the case of the first two partial derivatives in ([7]) are 

dtU = {U-krWkdti') + {U-kdtr){Uki') (19) 

and 

To recover an important property of the Wigner transform in the correspon- 
dence between the Liouville equation and TDSE, let us assume that V" = 0, 
so that kV in ([7]) can be written in the form 

kv'u = ^[{u-kr){UkVii^) - {u-kvr){Ukn ■ (21) 

Replacing ([19]), ([201) and (IH]) in ([7]) we obtain 

iU-kr){UkAiP) + {U^kA*r){Uki^) = , (22) 
where A is the linear operator 

A^d,-^dl + ^V . (23) 

Therefore, defined by (fT7|) is a solution of ([7]) if Aip = 0, or 

2 

iadtiP = HiP , H = -^dl + V , (24) 

formally the same as TDSE for the complex wave-function -0. A similar 
result can be obtained in the case of a charged particle in uniform magnetic 
field, presented in Appendix 2. 

For the harmonic oscillator potential V{x) = muj'^x'^/2, (j24p has both 
stationary (Hipn = Enipn) and nonstationary solutions, known as Glauber 
coherent states 

^c^^^^-^ ^ ^'^•^l^-iLvt/2^-a(x-u)'^/2+iav(x-u/2)/muj ^ ^25) 

vr 
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where a is a constant while u and v satisfy the equations of motion 

V 

dtu = — , dtv = —muju . (26) 
m 

By the Wigner transform (jlSp . '0^ yields the classical distribution 

/^c(x,p,t) = ^^e-"(^-")'-"(P-^)'/(™'^)' , (27) 

representing a Gaussian with a fixed width in the phase-space, oscillating 
along a classical trajectory. 

The limitation to polynomial potentials of degree at most 2 (square well, 
uniform field, harmonic oscillator), presumed above to derive (I24p . reflects 
the van Hove theorem on the validity of the canonical quantization [14] . 
In general, if V" ^ and Tp{x,t) is a solution of (f24l) . then f^{x,p,t) 
satisfies a modified Liouville equation [15]. For instance, if V{x) is a quartic 
polynomial, the modified equation has the form [L6l 

dtU + ^d^U-V'dpf^ = -^V"'d'pU ■ (28) 

One should note though that by considering f^{x, k,t) as a "matrix" element 
Pab = 'ip{xa)'ip*{xb) of the density operator p between Xa = x+ak/2 and x^ = 
X - ak/2, then for k^Q,V' = {V{xa) - V{xb))/{xa - x^), akV'f^ = [V, pU 
and ([7]) takes the form of the quantum Liouville equation iadtPab = [H, p]ab- 



3 The thermal sound 

The equations of motion for a nonrelativistic Brownian particle are 

dtx = ^ , dtp = -V' + m-l- (29) 
m m 

where 7 is the friction coefficient and 1^, —^p/m, denote the external force, 
respectively the backreaction produced by the thermal environment. 

In the bilinear coupling model |17J ^ is a function ^{£.Q,t) of time and 
the variables 8q describing the microscopic structure of the environment at 
the initial moment t = 0. If J-{£, T) is the distribution function of the 
environment in thermal equilibrium at the temperature T, then 

« m »= J dn£,T{£o,T)a£o,t) = (30) 
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and 

« m^{t') »= / dne,T{£o, T)^{£o, mSo, t') = kBTT{t - t') , (31) 

where T{t) is the memory function. The phase-space average << ... >> 
is presumed to be the same as the average < .. > over an ensemble of 
independent Brownian trajectories. Thus, the statistical significance of the 
distribution J-{£, T) appears indirectly, by the correlation function (j3ip . For 
the Brownian particle described by (j29p . this is specified by T{t) = 2-y6{t). 

Because of friction, a volume element = 6x5p shrinks according to 
dt5Q = —jSQ/m, and the particle number conservation dt{fSO,) = implies 
dtf = /m. The definition of the total derivative dtf 



d[_ 

dt 



dtf + idtx)d,f + {dtp)dpf , (32) 



and (j29p show that the local dynamics of / is expressed by 

9t/ + ^9,/-(y'-e + 7-)9p/ = ^/ , (33) 
m m m 

or 

dtf = --d.f + dp{V'-^ + j^)f . (34) 
m m 

Due to the noise, the statistical ensemble describing the particle is altered, 

so that the actual distribution function is the average < f > {x,p,t) of 

f{x,p,t) at each phase-space point {x,p). The transport equation satisfied 

by < / >, obtained from (j34p by average over an ensemble of solutions 

f{x,p,t), is the classical Fokker-Planck equation 

dt<f>+^d,<f>-dp{V' + j^)<f>=jkBTd^<f> . (35) 
m m ^ 

When describes an environmental particle, its outcome < / > should repro- 
duce self-consistently T{£,T). At zero temperature, (f35]) admits solutions 
of the form ([TT]) . with n{x,t) and S{x,t) provided by (fT3]) . respectively the 
modified Hamilton-Jacobi equation 

dtS + ^^^+V + ^S = , (36) 
Zm m 



associated to the classical equations of motion with linear friction jl8j . 
At finite temperature it is convenient to use the Fourier transform 

< f>ix,k,t) ^ j dp e*^^' < / > (x,p,t) , (37) 
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and the corresponding transport equation 

dt<f> --dlu <f> +k{iV' + ^dk) < f >= -jkBTk^ <f> . (38) 
m m 

A series expansion of < / > (x, A;, t) si k = 0, 

< / > (x, A;, t) = n(x, t) + imj{x, t)k 

- 2me{x, - i{2mfxix, + ... (39) 

provides the ensemble averages of probabihty ([8]), current Q and kinetic 
energy ([TO]) densities, denoted for simphcity in ([39]) also by n, j and e. 
When ()39p is replaced in (jSSp . the consecutive terms of the expansion yield 
the system of coupled equations 

dtn = (40) 
mdtj = -V'n - 7j - 2da:e (41) 

dte = ^^n-V'j-'^e-2d,X , - (42) 
m m 

If the cubic term x{x, t) = J dp < f > / (2m)^ in (j42]l is neglected, then 
thermal equilibrium {dte = 0) is reached when 

ksT mV . 
e= n 7 . (43) 

In this stationary state corresponds to dtj = and 

2 V V m 
j = —d^e n = -Dd^n n + -^d^{V'j) (44) 

^ ^ 7 7 

where D = ksT/j is the diffusion coefficienlH. At strong friction the last 
term in ()44|) can be neglected, and ()40|) becomes the Smoluchowski equation 

dtn = d^{Dd,,n+—n) . (45) 
7 

In states of nonequilibrium at low temperature, the kinetic equation ([5]) 
provides the velocity of zero sound [19j. Similarly, applied to a nonstationary 
state, (go]) and dH]) yield 

d?n = -dltj = -dle + -d,{V'n)-^dtn . (46) 
mm m 



^If dxt = too, then j = —V'n/j and according to (|43|) . e is the sum between a drift 
component mj^ /2n and a thermal component nkBT/2. 
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If y = and e is given by (|43p . then (j46p becomes a wave equation for the 
probabihty density, 

a^n = w^a^n - —dtn , (47) 



where Vs = ^Jk^Tjm is the sound velocity provided by Newton's formula. 

It is interesting to remark that if of (|17p is expanded in the form (|39p , 
then i = ndxS/m. e = mj^/2n + a'^[{dx\'^)^ — ^Jnd1.^/ri\l Am^ so that when 
T = 0, 7 = 0, (HDD and dni) yield 



A„ = -8.i. „S.|8.S + %3^-|lS^ + V] = 0. (48, 

close to ([13]), ([M]) and equivalent to the TDSE jMl). Moreover, if < > 
(x, fc, t) is expressed as a matrix element < p >ab (t) of the average density 
operator < p > between Xa = x + ak/2 and Xb = x — ak/2, then for A; 
()38p takes the form of the quantum Fokker-Planck equation [20j . 



4 Summary and conclusions 

Probability waves in the configuration space are related to certain solutions 
of the Liouville equation which are coherent, in the sense that they keep in 
time the same functional dependence on the coordinate and momentum. Of 
course, this kinematical property alone does not explain how such distribu- 
tions are formed, but only shows that once created, they are stable. 

Two generic examples of classical, functional coherent states, are pre- 
sented in Sect. 2. The first example concerns the action distributions (jlip. 
localized in the momentum space at the value given by the generating func- 
tion S. For the related action wave the current (j) and kinetic energy (e) 
densities are expressed in terms of only two functions, the probability den- 
sity n and S, provided by the continuity, respectively the Hamilton-Jacobi 
equations. These equations can also be obtained from a least-action princi- 
ple in which n and S are canonically conjugate (Appendix 1). 

The second example is derived from the first, presuming the existence 
of a minimum distance in the coordinate space which decreases as a func- 
tion of the inertial parameter. Worth noting, this elementary distance may 
not necessarily result from a lattice structure of the phase-space, but take 
into account a minimum interval of time [5]. It is shown that if the space 
derivative of S in the action distribution is written as finite difference, then 
one obtains a Wigner-type distribution psp . expressed in terms of a single 
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complex functional ip of n. However, by contrast to the action distribu- 
tions, the coherence of the Wigner distributions is maintained only by a 
harmonic oscillator, uniform field, or square well potential. For the har- 
monic oscillator, beside the static probability distributions can be found 
Gaussian wave-packets in phase-space (f27|) . oscillating along the classical 
orbit. Such solutions could be useful for instance to describe a macroscopic 
Bose-Einstein condensate in a harmonic trap [21]. 

The evolution of the average phase-space distribution function for a 
Brownian particle is discussed in Sect. 3. Containing stochastic and dis- 
sipative forces, the equations of motion ([29|) yield the modified Liouville 
equation p4p and then, by ensemble average, the Fokker-Planck equation 
(j35p . The Fourier transform in momentum of the average distribution func- 
tion provides the localization probability (n), current (j), and kinetic energy 
(e) densities, while the pure distributions discussed in Sect. 2 may still be 
useful to obtain solutions in the limit of zero temperature. However, in 
general (jSSp describes the irreversible increase of the entropy, and the func- 
tional relationship between n, j and e is specified only by the system (|10|) . 
(I4ip . (I42p . In free space at thermal equilibrium this system reduces to the 
sound waves equation (j47p . while for a Wigner distribution at 7,T zero it 
provides the Schrodinger equation (H8]l . As both equations arise in the same 
framework, the result of thermal averaging in nonrelativistic quantum sys- 
tems could be a specific form of quasiclassical behaviour, expressed by a 
transition from complex (tp) to real (n) probability waves, as an intermedi- 
ate stage between decoherence and complete dissipation. One should note 
though that at strong friction this stage can be suppressed, and in the pre- 
vious calculations on atomic tunneling at finite temperature [6], it was not 
observed. 

The assumption << >>~ 6{t — t') in ([3T|) . of 5-correlated noise, 

is reflected by the second-order derivative of the distribution function in the 
Fokker-Planck equation (j35p . This term is essential to retrieve the nonrel- 
ativistic wave equation (I47p . but is not suitable to simulate effects of the 
quantum fluctuations, such as the third-order derivative dpf^ in (j28p . The 
emergence of the quantum dynamics out of a discontinuous character of the 
inertial motion, rather than by the action of random forces, might provide 
an alternative worth of further consideration. 
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5 Appendix 1: Action functional for action waves 



The coupled equations p3|) , p4|) describing action waves can also be derived 
using the variational principle Sn,s-^ = for the action functional ^ 



A[n, S] = - I dt [ dx n\dtS + -^^^ + 1^1 (49) 
J J 2m 

with respect to the local variations of the fields n(x, t) and S{x^ t). The vari- 
ational formulation shows that the coupled continuity and Hamilton-Jacobi 
equations describe an infinite-dimensional Hamiltonian system, having the 
generating function S{x, t) and the density n{x, t) as conjugate variables. 
The symplectic form and the Hamilton function of this system are, respec- 
tively 

uj= [ dx {dnA dS) (50) 



and 

HA = ldx + V] . (51) 



In this framework, the conservation of the particle number expressed by (jl3p 
reflects the " gauge-invariance" of the variational equation Sn,s^^ = with 
respect to the change of S by adding an arbitrary function of time. 

It is interesting to note that (j49p . (jSOp and (jSip define a Hamiltonian 
system in the Hilbert space TC generated by = \/nex.p{iS/a), while the 
equality 

dx {dn A dS) = -ia / dx {d^* A dip) (52) 



proves that iv coincides with the symplectic form |22j induced by the complex 
structure of TC. 



6 Appendix 2: Quantum waves in magnetic field 

In the three-dimensional space, an electric charge in the magnetic field B = 
V X A, is described hy H = {p - qA)^/2m + V, and ([5]), dH) take the form 

dtf+^^-^-Vf-[V-Vp ,y]/+ ^^P~^^^ [V-Vp ,A]/ = , (53) 
m m 

respectively 

2 

9t/--V-Vfc/+i[k-V,y+^A2]/--(A-V+V,.-[k-V,A])/ = , (54) 
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where V = d/dv, Vp = d/d-p and = d/d\s.. 

For a quantum nonrelativistic scalar particle the coherent solutions of 
(j54p have the form 

/^(r,k,t)^V'*(r-— ,t)V'(r + — ,t) = (t/-V)(C/^) (55) 

with (j = e^'^'^/^. Presuming that A = (B x r)/2, with B constant, and 
(k • V)'^y = 0, then for the first three terms of ()54p we obtain 

2 

St/v, - -V • Vkh + ^[k • V, y + :^A2]/^ = (56) 



where 



(f7^V*)(f^AiV^) + {U-^Ali^*)iUi') 

Ai^d,-'^V' + Uv + ^A') . (57) 
Zm n zm 



Because V • A = 0, Vfe • [k • V, A] = [k • V, A] • Vfe, and in the last term 
of (ISl 

(A-V + [k- V,A]- Vfc){7V' = f7(A-V)V' • (58) 
Thus, with dMl) and (IMI), (IMl) takes the form 

{u-^r){uh^i') + {u-^kir){u^) = ^ , (59) 

where A2 = Ai — qA • V/m. This equation reduces to A2V' = 0, or ihdtip = 
Hip, with 

zm 

For a wave-function ^ with two components, describing a particle with 
spin 1/2, ()54p may have coherent solutions of the form 

U{r,k,t) ^ {U-'^)HUs^) (61) 

where f/^ = e^^'''^^^^'®)/^, and s is the spin operator. With Ug instead of U, 
to the first order in k ()56p remains the same, but (j58p is changed to 



(A•V + [k•V,A]•Vfc)^7,^' = ^7,(A•V + ^s•B)^ , (62) 



so that (1601) becomes 



Hs = H- — s • B . (63) 
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Therefore, the solution ^'(r, t) of the TDSE ifidt^ = Hg^ for a particle with 
spin corresponds to a coherent solution 

/*(r,P,t) = ^ j d'k e-*'^ P(?7-i^)t([7,M/) (64) 

of the classical Liouville equation (I53p . One should note though that in ()63p 
the magnetic moment q7i/4m is half the value provided by the relativistic 
Dirac equation, and a complete description of a charged particle with spin 
should include the internal degrees of freedom also in the Liouville equation. 
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